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While microscopic laws of physics are invariant under the reversal 
of the arrow of time, the transport of energy and information in 
most devices is an irreversible process. It is this irreversibility that 
leads to intrinsic dissipations in electronic devices and limits the 
possibility of quantum computation. We theoretically predict that 
the electric field can induce a substantial amount of dissipationless 
quantum spin current at room temperature, in hole doped semicon- 
ductors such as Si, Ge and GaAs. Based on a generalization of the 
quantum Hall effect, the predicted effect leads to efficient spin injec- 
tion without the need for metallic ferromagnets. Principles found in 
this work could enable quantum spintronic devices with integrated 
information processing and storage units, operating with low power 
consumption and performing reversible quantum computation. 

Our work is driven by the confluence of the important technological goals of quantum 
spintronics (0,0) with the quest of generalizing the quantum Hall effect (QHE) to higher 
dimensions. The QHE is a manifestation of quantum mechanics observable at macroscopic 
scales. In contrast to the most transport coefficients in solid state systems, which are 
determined by the elastic and the inelastic scattering rates, the Hall conductance on in 
QHE is quantized and completely independent of any scattering rates in the system, where 
the transport equation is given by j a = autapEp ( j a and Ep (a, (3 = 1, 2) are the charge 
current and the electric fields respectively, e a p is the fully antisymmetric tensor in two 
dimensions ). While dissipative transport coefficients are expressed in terms of states in 
the vicinity of the Fermi level, the non-dissipative quantum Hall conductance is expressed 
in terms of equilibrium response of all states below the Fermi level. The topological 
origin of the QHE is revealed through the fact that the Hall conductance can also be 
expressed as the first Chern number of a U(l) gauge connection defined in momentum 
space (0)- Recently, the QHE has been generalized to four spatial dimensions (4)- In that 
case, an electric field E v induces an SU{2) spin current f = 1,2,3,4, i = 1,2,3) 

through the non-dissipative transport equation f = ar/^ u E u , where rj 1 is the t'Hooft 
tensor, explicitly given by 77* = e i/LU ,4 + 5^5^ — 5 iv 5^ and a is a dissipationless transport 
coefficient. The quantum Hall response in that system is physically realized through 
the spin-orbit coupling in a time-reversal symmetric system. At the boundary of this 
four dimensional quantum liquid, when both the electric field and the spin current are 
restricted to the three dimensional sub-space, the dissipationless response is given by 

j) = a s e^ k E k (1) 

This fundamental response equation shows that it is possible to induce a purely topological 
and dissipationless spin current by an electric field in the physical, three dimensional space. 

We consider a realization of this electric field induced topological spin current in 
conventional hole-doped semiconductors. In a large class of semiconductors, including Si, 
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Ge, GaAs and InSb, the valence bands are four-fold degenerate at the T-point (see e.g. 
Fig. 1). The effective Luttinger Hamiltonian (J3J) for holes is given by 

^o = |^((7i + ^7 2 )fc 2 - 2 72 (k-S) 2 ), (2) 

where Si is the spin-3/2 matrix. We take the hole picture, and reverse the sign of the 
energy. Good quantum numbers for this Hamiltonian are the helicity A = /i _1 k- S/k, and 
the total angular momentum J = hx. x k + S. This kinetic Hamiltonian is diagonalized 
in the basis where the helicity operator A is diagonal, and the eigenvalue is given by 
e\ (k) = ^-(71 + (§ — 2A 2 )7 2 ) = For a given wave vector k, the Hamiltonian © 

has two eigenvalues, e H (k) = e A=±3/2 (k) = :a ^ 11 h 2 k 2 = |^ and e L (k) = e A=± i /2 (k) = 

11 2m 2 = forming Kramers doublets. They are referred to as the light-hole (LH) 
and heavy-hole (HH) bands. In semiconductors with zincblende structure, such as GaAs, 
inversion symmetry breaking causes an additional tiny splitting in the LH and HH bands. 
We can neglect it when the temperature is much higher than this splitting. The band 
structure of semiconductors deviates from the spherical to the cubic symmetry. We also 
neglect this effect for simplicity, because physics described below are not so much affected 
by it. 

We shall consider the effect of a uniform electric field E. Our full Hamiltonian is thus 
given by H = H + V(x), where V(x) = eE ■ x, and — e is the charge of an electron. 
We assume that the split-off band is totally occupied. We first define a 4 x 4 unitary 
matrix U(k) which diagonalizes the kinetic Hamiltonian H Q . U(k) is defined by £/(k)(k • 
S)L/t(k) = kS z . In the spherical coordinates where k = k(sin8 cos0, sin# sin 0, cos#), 
U (k) can be expressed as U (k) = exp(i9S y ) exp(i<pS z ). Under this unitary transformation, 
the new Hamiltonian H = U(k)HU'(k) becomes 

H = ^ (71 + \l2 - 2 l2 Sl) + (k)V(x)tft(k) (3) 

Eigenvalues of S z physically describe the helicity A = ftT k ■ S/k in the original basis. 
The kinetic part Hq now becomes diagonal, in the representation where S z is diagonal. 
Because x = id^, the potential term becomes V(D), where the covariant derivative D is 
defined by D = id^ — A and A = —ill (k)<9 k £/"l"(k). As A is a pure gauge potential, there 
is no curvature associated with it. Up to this point, the transformation is exact. We 
now consider adiabatic transport and make a corresponding approximation. As is usually 
assumed in the transport theory, we neglect the interband transitions, i.e. the off-block- 
diagonal matrix elements of A connecting the LH and HH bands. Then we arrive at a 
non-trivial adiabatic gauge connection A (online supporting text), which takes a block- 
diagonal form in the LH and HH subspace. As each band is two- fold degenerate, the gauge 
connection is in general non-Abelian. However, A has no matrix elements connecting the 
A = 3/2 and A = —3/2 states in the HH band, because the gauge field A only connects 
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states with helicity difference AA = 0, ±1. Therefore, the non-Abelian structure is only 
present in the LH band. For simplicity of presentation, we shall first make an additional, 
Abelian approximation (AA), in which only the diagonal components in A are retained. 
Afterwards, we shall give our final results including fully the non-Abelian corrections. 

Within the AA, A is a diagonal 4x4 matrix in the helicity basis. As a band-touching 
point acts as a Dirac magnetic monopole in momentum space (0), each diagonal compo- 
nent of A(k) is given by that of a Dirac monopole at k = 0, with the monopole strength 
eg given by A. The associated field strength is given by 

F ij = i[D i ,D j ] = e ijk X^ (4) 

The effective Hamiltonian takes the form 

= + V,x). ,5) 

Henceforth, x^ denotes a covariant derivative in momentum space: Xi = -D, = id/dki — 
Aj(k). Note that the definition of x, has changed by projecting the original Hamiltonian 
H onto the HH or LH band. While H eS seems to be trivial, its non-trivial dynamics is 
revealed through the non-trivial commutation relations 

[ki, kj] = , [Xi, kj\ = i5ij , = —iFij. (6) 

Such situation also happens in the Gutzwiller projection of the SO (5) model Q. It 
also resembles the non-trivial commutation relation between the position operators of 
a two-dimensional-electron-gas projected onto the lowest-Landau- level (0), where Fjj = 
Btij, and B is the external magnetic field. This general algebraic structure, called "non- 
commutative geometry", also underlies the four-dimensional QHE model @. In our 
present context, the non-commutativity between the three-dimensional coordinates arises 
from the magnetic monopole in momentum space, and it is a natural generalization of the 
QHE to three dimensions. 

The equation of motion for holes can be derived easily from Eqs. and El as 

hk- 

Tiki = eE u Xi = — - + Fijkj. (7) 
m x 

The last term, proportional to Fy, is a topological term, describing the effect of the 
magnetic monopole on the orbital motion. It represents a "Lorentz force" in momentum 
space, making the hole velocity non-collinear with its momentum, in contrast to the usual 
situations. In fact, if we interchange the roles of x and k in this term, it becomes the 
Lorentz force for a charged particle moving in the presence of a magnetic monopole in real 
space. This set of equations can be integrated analytically (supporting online text), and 
the resulting trajectory is shown in Fig. 2. The hole motion in real space obtains a shift 
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perpendicular to S. This shift is analogous to the deflection of a charged particle by a 
magnetic monopole, in a direction perpendicular to the plane spanned by its position and 
velocity vectors (0). It causes a spin current perpendicular to both E and S. For example, 
for E parallel to the +z direction, the spin current for each band at zero temperature, 
with spin parallel to the x axis, flowing to the y direction is given by 

pF 

^ = 36^ (9 ^ + ^' (8) 

which is obtained by summing contributions from all the filled states. Here we assumed 
that the equilibrium momentum distribution is attained by the random impurity scat- 
tering which causes the charge relaxation. Note that Eq. [7| describes only the ballistic 
motion and scattering by random impurities would lead to additional contributions to the 
spin current. As one can see from the detailed discussions in the supporting online text, 
these extrinsic effects are not only small, but also scale with a higher power of kp ~ n 1//3 , 
where n is the hole density. Therefore, by plotting Us/n 1 / 3 against n, and extrapolating 
to the limit of n — > 0, the constant intercept would uniquely determine our predicted 
dissipationless spin conductivity. 

It is worth noting that this AA becomes exact in zero-gap semiconductors, e.g. a-Sn. 
In this class of materials, the bottom of the conduction band and the top of the valence 
band correspond to the LH and HH bands in other semiconductors like GaAs. These two 
bands touch at k = 0. In this case, p-doping introduces holes only into the HH band, and 
the A A becomes exact. 

The electric-field-induced spin current can also be understood in terms of the conser- 
vation of the total angular momentum J = fox x k + S. As remarked earlier, J commutes 
with Hq. When E is parallel to the z-direction, J z also commutes with the potential. 
Therefore, substituting S = Afok = \hk/k, we obtain 

j z = fo(x x k) 2 + fo(x x k) 2 + Afok 2 = (9) 

The second term, representing the torque, vanishes in our case since k points along the 
z direction. The first term fo(x x k) z vanishes in usual problems; however, it does not in 
our case, due to the non-collinearity of the velocity and the momentum. Furthermore, the 
first term, describing the time derivative of the orbital angular momentum L = fox x k, 

is proportional to the spin current. The third term Afok 2 describing the time derivative 
of the spin angular momentum S, can be easily evaluated from the acceleration equation 
in Eq. [7| Therefore, we see that the conservation of the total angular momentum Eq. |U] 
directly implies the spin current Eq. |HJ The spin current flows in such a way that the 
change of L exactly cancels the change of S. 

We now discuss the correction due to the non-Abelian nature of the gauge connection 
of the LH band. Remarkably, even though the gauge connection is non-Abelian, the 
associated field strength is Abelian, and gives a correction factor of (—3), compared 
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to the AA (fT^fTU fT ^) . The equation of motion modified accordingly agrees with that 
obtained by generalizing the wave packet formalism to the non-Abelian case. This 
non-Abelian correction gives the following result for the spin current: 

fy = ^ k *-k h F ) = ^ s E z . (10) 

Here we defined a s to have the same dimension as the electrical conductivity, to facilitate 
comparison. The spin current equation is rotationally invariant, with the covariant form 
given in Eq. ^ and is the central result of our paper. In contrast with similar effects 
QUI, this spin current has a topological character; the spin conductivity a s in Eq. ^ 
is independent of the mean free path and relaxational rates, and all states below the 
Fermi energy contribute to the spin current, where each contribution is purely determined 
by the gauge curvature in momentum space, similar to the QHE (3). Assuming the 
hole density n = 10 19 cm~ 3 , the mobility of the holes at room temperature in GaAs is 
/i = 50cm 2 /Vsec (f73J). and the conductivity is a = en\i = 80fi _1 cm" 1 . On the other 
hand, the spin Hall conductivity a s in Eq. ^TO] is estimated as <j s ~ 80fi _1 cm _1 , being of 
the same order with a . For lower carrier concentration, a s becomes larger than a; for 
n = 10 16 cm -3 , we have o = 0.6f2 _1 cm _1 and a s = 7f2 _1 cm _1 . At finite temperature, 
Eq. El is modified only through the Fermi distribution function n x (k). Since the typical 
energy difference between the LH and HH bands at the same wavenumber is about O.leV, 
which largely exceeds the energy scale of the room temperature ~ 0.025eV, our predicted 
effect remains of the same order even at room temperature. 

We remark that the non-dissipative spin transport equation Eq. ^does not violate the 
time-reversal symmetry T. Our microscopic Hamiltonian H, the electric field E and the 
spin current are all T invariant. Therefore, the electric field and the spin current can be 
related by a T symmetric, dissipationless transport coefficient a s . This situation is to be 
contrasted with the Ohm's law. As the charge current is odd under T, while the electric 
field is even, they can only be related by a T antisymmetric, dissipative transport coeffi- 
cient, namely the charge conductivity. One of the main objective in quantum computing 
is to achieve reversible computation fliTpiffi . From the above analysis, we see that there is 
a fundamental difference between the ordinary irreversible electronics computation based 
on the Ohm's law, and the reversible spintronics computation based on Eq. ^ The time 
reversal symmetry property encoded in Eq. ^ could provide a fundamental principle for 
the reversible quantum computation. 

This spin current is also useful for spin injection into semiconductors. While effective 
spin injection is necessary for spintronic devices, it has been an elusive issue (pj. Usage of 
ferromagnetic metals is not practical because most of the spin polarizations are lost at the 
interface due to conductivity mismatch between metal and semiconductor (fT^ I ^) . Spin 
injection from ferromagnetic semiconductors such as Gai_ x Mn x As has been successful $21\ 
HUH3J). Nevertheless, T c is at most 110K for Gai_ x Mn x As, still too low for practical use at 
room temperature. Thus it is desirable to find an effective method for spin injection. The 
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electric-field-induced spin current serves as a spin injector, because it creates a spin current 
inside the semiconductor. One might worry that the short relaxation time r s ~ lOOfsec 
(24) of the hole spins. This shortness of t s is because the strong spin-orbit interaction in 
the valence band combines the relaxation of momenta and spins (JHJ). Actually most of 
the efforts on the spintronics in GaAs have been focused on electrons in the conduction 
bands, having much longer spin-relaxation time (~ lOOpsec) (|i?6j). Nevertheless, our 
spin current is free from such rapid relaxation of spins of holes, because it is a purely 
quantum mechamical effect with equilibrium spin/momentum distribution. Only when 
the spin/momentum distribution deviates from equilibrium, e.g., in spin accumulation at 
boundaries of the sample, the rapid relaxation of hole spins becomes effective. 

One can consider following experimetal setups for detection of the constant spin sup- 
ply from p-GaAs. When the electric field is applied along the ^-direction and the electric 
current J z is induced, the s x -spin current j® will flow along the y direction. One possi- 
bility is to see the spin-dependent electric transport through a ferromagnetic electrode 
with the magnetization M along indirection attached to the positive-?/ side of the sam- 
ple. With a lead connecting this electrode and the other (negative-?/) side of p-GaAs as 
shown in Fig. 3 A, one should see a change of the electric current / depending on the 
direction of M. The ratio of / when M is along indirection, I(+x)/I(— x), is expected 
to be well larger than unity. For the ferromagnetic electrode, ferromagnetic metals are 
not efficient, because of the conductance mismatch <\2CQ . Instead, ferromagnetic semicon- 
ductors will be suitable. Another possibility is to measure circular polarization of light 
emitted via recombination with electrons. This can be achieved by a similar experimental 
setup in Ref. (p?Tj). replacing (Ga,Mn)As by p-GaAs, where the quantum well structure 
of (In,Ga)As is sandwiched by p-GaAs and n-GaAs (Fig. 3B). The spin current injected 
along the ^/-direction will be recombined with the electrons supplied from the n-GaAs in 
the (In,Ga)As quantum well. 

When the system is not connected to the leads along the y-direction, spins accumulate 
near the edges of the sample. This spin polarization can in principle be measured by the 
Kerr rotation. The spin distribution is determined by a balance between the spin current 
supply and the spin relaxation. At room temperature, because r s = lOOfsec {24) is rather 
short, the area density of spin accumulation at the sample surface, jyT s , is too small 
to be observed. However, there are several ways to make t s longer. One is to lower 
the temperature. Another is to inject the spin into n-GaAs through the p-n junction as 
demonstrated recently ( pffiUffi) . The spin lifetime of electrons in GaAs is lOOpsec, 10 3 
times longer than that of holes. Thereby the spin current can be detected by the Kerr 
rotation through surface reflection. 
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Figure 1: Approximate band structure of GaAs. We neglect the small splitting due 
to inversion symmetry breaking. We also neglect the anisotropy of the bands. The 
conduction band (CB) is two-fold degenerate. The valence band consists of the heavy 
hole (HH), the light hole (LH), and the split-off (SO) bands, each of which is two-fold 
degenerate. We consider the p-GaAs, and the Fermi momentum for each band is labelled 
as kp and k F respectively. The Fermi energy shown in the figure corresponds to n — 
10 19 cm~ 3 . 
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Figure 2: The real-space trajectory of the hole obtained by solving Eq.(J7J), The electric 
field E = e _1 VV is parallel to the +z direction. Due to the non-commutative relation 
between the components of the position operator x% y i.e., = —iFij with Fy being 

the gauge curvature defined in the momentum space, the hole obtains a transverse velocity 
whose direction depends on the helicity A = k • S/k indicated by the thick arrow. The 
broken line is parallel to k. 




Figure 3: Experimental setups for the detection of spin current induced by an electric 
field. (A) Detection by attaching a ferromagnetic electrode. The dependence of current 
/ flowing into the electrode on the direction of the magnetization M is to be measured. 
(B) Detection by measuring the polarization of the emitted light in the quantum well of 
(In,Ga)As. Right and left circular polarization will be switched when the direction of the 
external electric field is reversed. 
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Luttinger Hamiltonian for semiconductors 

Here we briefly summarize basic aspects of the Luttinger Hamiltonian for semiconductors. 
In isolated atoms, spin-orbit coupling leads to a splitting of the p orbitals into four-fold 
degenerate P3/2 and two-fold degenerate P\i2 levels. In a large class of semiconductors, 
including Si, Ge, GaAs and InSb, the P3/2 levels form the top of the valence bands, 
which are separated from S-like conduction bands. Therefore, the valence bands are four- 
fold degenerate and the conduction bands are two-fold degenerate at the T-point (see 
e.g. Fig. 1). The form of the effective Hamiltonian including spin-orbit coupling can 
be determined from symmetry arguments alone. In crystals with inversion symmetry, 
band energies depend quadratically on the momentum away from the center; therefore, a 
rotationally invariant Hamiltonian can only contain two possible terms, k 2 and (k • S) 2 , 
where S is the 2x2 Pauli matrix for the two-fold conduction band and the 4x4 spin matrix 
for the four-fold valence band. As a square of each Pauli matrix is the identity matrix, 
to the lowest order, there is no spin-orbit coupling in the conduction band with inversion 
symmetry. In the presence of inversion symmetry breaking, either due to intrinsic crystal 
symmetry or in hetero-junction systems, the Dresselhaus and the Rashba terms for spin- 
orbit coupling are possible (S1-S3). On the other hand, in the four-fold valence band, 
spin-orbit coupling affects the band structure, even without inversion symmetry. As the 
square of each component Si {i = 1,2,3) for S = 3/2 is non-trivial, these two terms 
combine to form the effective Luttinger Hamiltonian (S4) for holes: 

Ho = ^ ((7i + \l2)k 2 - 2 72 (k ■ S) 2 ) , (SI) 

where we take the hole picture, and reverse the sign of the energy. The eigenvalues are 
classified by the helicity A = fc _1 k • S/k; A = ±1/2 corresponds to the LH band, and 
A = ±3/2 to the HH band. 



Details for calculation of A(k) 

The spin-3/2 matrices Si are written as 

\ / 

1 

1 



S x 



( yj 

2 

2 



V 



2 



yj 
2 



S,, 



2 1 



V 



2 1 



2 1 



j S z 



(S2) 
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By using U(k) = exp(i9S y ) exp(i(j)S z ), the gauge field A = — iU(k)c\W(k) is calculated 

as 



A-dk 



■| cos Odtp 



?§(sm6d(p-id6) 



^(sm9d<p + id9) 
— \ cos 9d(f 
sin 9dip — id9 



\ 



sin 9d(p + id9 
\ cos 9d(p 
^(sm9dip-id9) 



^-(sm9dip + id9) 

\ ^ (sin 9d(p - %d9) § cos 9d(p J 

' " (S3) 

where the first and fourth columns correspond to the HH, while the second and third 
ones correspond to the LH. As A is a pure gauge potential, its curvature Fij = i[Di,Dj] 
vanishes. In considering adiabatic transport, we neglect the interband transitions, i.e. the 
off-block-diagonal matrix elements of A between the LH and HH bands. Then we arrive 
at a non-trivial adiabatic gauge connection 



A'-dk 



/ — | cos 9d(p 



V 



— \ cos 9 dip 
sin 9dip — id9 



sin 9dip + id9 
\ cos 9dip 



\ 



| cos 9dip J 



(S4) 



which takes a block-diagonal form in the LH and HH subspace. As the states in each 
band are two-fold degenerate, the gauge connection is in general non-Abelian. However, 
the non-Abelian structure is only present in the LH band; namely, in the LH subspace the 
matrix A' has off-diagonal components, while in the HH subspace A' is diagonal. This 
is because the gauge field A only connects states with helicity difference A A = 0,±1, 
and does not connect A = ±3/2 states in the HH band. By employing the Abelian 
approximation (AA) we neglect the off-diagonal elements of A', and we get thereby 



A' Abelian = -S 2 COS 6>V0, 



(S5) 



which is diagonal in the basis where S z is diagonal (S5). As has been recognized in 
(S6), the degeneracy point acts as the source of the gauge field, i.e. Dirac magnetic 
monopole. The present form of A' corresponds to a Dirac monopole at the origin of the 
momentum space, with the monopole strength eg given by S z , in accordance with (S6). 
To return to the original basis where the helicity A becomes diagonal, we make the gauge 
transformation 



A 



Abelian 



» = C/t( k )A' Abelian (k)C/(k) +;^M[/(k) 



dk 



(S6) 



The monopole strength is then expressed as eg 
strength is given by 

=i[Di, Dj] = € ijk \ 



A. The associated magnetic field 

h 



k 3 



(ST) 
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We now discuss the correction due to the non-Abelian nature of the gauge connection 
of the LH band. Remarkably, even though the gauge connection is non-Abelian, the 
associated field strength is Abelian, and is given by (S5) 

^ = ^ (2A 2 - I) (S8) 

This has a correction factor of (—3) only in the LH compared to the AA. As the spin 
current depends only on the field strength in momentum space, the spin current acquires 
an extra factor of (—3) only in the LH, compared with that within the AA. 



Real-space trajectory of the hole motion 

The equation of motion (Eq. 7) within the AA can be integrated analytically. When E is 
parallel to the +z direction, we get 

k x (t) = Ko, k y (t) = k y0 , k z (t) = k z0 + eE z t/h, 

u\ 1 hk z0 eE z 2 

z(t) =z + 1 + - — t , 

m\ 2m x 

hk x0 \k y0 eE z t/h + k z0 
x(t)=x + —t- y 



mx k 2 x0 + k 2 Jk 2 xQ + k\ ] + (eE z t/h + k z0 ) 2 



Hkyo \k x0 eE z t/h + k z0 

y(t) = y + —tt-t + 



mx k 2 x0 + k 2 y0 ^ k 2 Q + k 2 Q + (e Ez t/h + k z0 ) 2 ' 

The last terms in x(t) and y(t) represent a shift of the particle position, in a direction 
perpendicular to k, as presented in Fig. 2. As the directions of the spin is parallel 
(antiparallel) to k for A > (A < 0), the hole motion in real space obtains a shift 
perpendicular to S. This shift causes a spin current perpendicular to both E and S. 



Calculation of spin current 

The motion along the ^-direction in the equation of motion (Eq. 7) is free acceleration by 
the electric field. In reality, because this acceleration is suppressed by random scattering 
we should take into account the charge relaxation by random scattering as is discussed in 
details in the next section. For example, the spin current for each band at zero tempera- 
ture, with spin parallel to the x axis, flowing to the y direction is given by 

h -sjj/r^-T^nr-i^' (S9) 

•xL _ ^ „;^ x m x (ir\ - e tp V" nL ( k )^ _ eE z kp 

A=± 5 ,k K 
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where n A (k) is a filling of holes in the band with helicity A; n ±3//2 (k) = n H (k) and 
n ±1 / 2 (k) = n L (k) are fillings for the HH and LH bands, respectively. Charge relaxation 
is already included in this result by putting the Fermi surfaces as in their equilibrium 
position |k| = hp- We have taken into account that the expectation value of the electron 
spin s is one-third of that of S: s = |S. This is because the spin-3/2 matrix S is a 
summation of the spin angular momentum s with spin one-half and the atomic orbital 
angular momentum 1 with spin one (S4)- 

By including the non-Abelian correction in the LH band (Eq. IS8|) . the spin current is 
given by 

fy = ^tr^rWS x n\k) = ^tr^F yz (k)S x n\k) 




(Sll) 



Interestingly, in the Ioffe-Regel limit of kpl ~ 1 the size of the spin conductivity is 
comparable to that of the charge conductivity 

The formula, Eq. IS11[ is at zero temperature; for finite temperature, it is modified 
only through the Fermi distribution function n A (k). For example, at room temperature 
and n = 10 19 cm -3 , the nominal value of the energy difference between the LH and HH 
bands at the same wavenumber is O.leV, and it largely exceeds both the temperature 
(~ 0.025eV) and the energy scale for momentum relaxation h/r p ~ 0.006eV, where we 
estimated the momentum relaxation time r p to be of the same order (S7) as the spin 
relaxation time for holes r s ~ lOOfsec (S8). Thus the value of the spin current remains 
of the same order as in the zero temperature. 



Impurity scattering and relaxation 

Our equation of motion Eq. 7 describes only the ballistic motion, and the relaxation 
process due to the random scattering by impurities and phonons are not explicitly taken 
into account. The relaxation is essential to attain the thermal equilibrium, which we 
assumed to obtain Eq. 10. Our result is similar to the Karplus-Luttinger (KL) term (S9) 
in the context of the anomalous Hall effect (AHE). The effects of the random scatterings 
have been discussed in the context of AHE by many authors (S10-S14)', here we shall use 
some of these results to estimate the impurity contribution to the spin current. Kohn- 
Luttinger (Sll) and Luttinger (SI 2) derived the formula for AHE based on the rigorous 
treatment of the density matrix in the expansion of the impurity scattering strength v, 
reproducing all the contributions to a xy up the v order. The skew scattering contribution 
proposed by Smit (S10) is of the order of v" 1 , and hence is expected to be important in 
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the weak disorder case. However it is proportional to (V^ 3 ),; (V is the impurity potential, 
and () c is the cumulant average) and hence vanishes for the Gaussian distribution of the 
disorder potential. Furthermore, this skew scattering mechanism predicts pu oc p (p: Hall 
resistibity, p: diagonal resistivity), which is not usually observed experimentally including 
the doped GaAs. The contribution of the order of v contains two terms, both of which are 
independent of the scattering. One is the KL term, which is similar to our dissipationless 
spin current, and the other term is later interpreted by Berger (SI 3) as the side-jump 
contribution accompanied by the impurity scattering. In the following, we shall estimate 
their relative magnitudes and point out the important fact that they scale differently with 
the hole density n, and can therefore be uniquely distinguished in the low density limit. 

In the present case of GaAs, there are several evidences that the intrinsic AHE (KL 
term) is the dominant contribution. Jungwirth et al. (SI 5) showed that KL term ex- 
plains the experimentally observed a xy quantitatively for ferromagnetic (Ga,Mn)As and 
(In,Mn)As. Especially the difference between (Ga,Mn)As and (In,Mn)As is successfully 
attributed to the different ratio m^uj^HH of the effective masses of light- and heavy hole 
bands. Also one can estimate the contribution of side jump and skew scattering as follows. 
By modifying the expression for the AHE due to side jump mechanism (S13,S14), the 
side-jump contribution to the spin Hall conductance is given by 

2 \ 2 

^side jump = -^|A n (S12) 

where A c = h/mc = 3.86 x 10~ n cm is the Compton wave length of free electron, and a 
is the enhancement factor for the Bloch electrons which is estimated to be around 10 4 . 
In terms of this estimation one can estimate its ratio to the spin Hall conductance a s in 
Eq. 10 as 



0-side jump 

= UL\A C KF^ 2 

(7, 

120 



a(\ c k F f (S13) 



which is around 10 even for the density as high as n = 10 cm . Therefore it is 
expected to be negligible. As for the skew scattering contribution, it is at most of the 
order of 



_skcw scattering 
° s 



a(X c k F ) 2 f-^ (S14) 



which is again much smaller than unity for the reasonable value of y ^f 1 ~ 10. Therefore 
we have the good reason to expect that our <j s in Eq. 10 gives the dominant contribution, 
and hence the spin current estimated in this paper can be observed in hole doped GaAs. 

From Eqs. IS12US13l and lS14l we see that the side jump and the skew scattering contri- 
butions scale with the hole density like n and n 4//3 respectively, while our dissipationless 
spin conductivity a s scales like n 1 / 3 . Therefore, one can plot the experimentally observed 
total spin conductivity as crj/n 1//3 versus n, and extrapolate to the limit of n — > 0. Only 
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the dissipationless spin conductivity contributes to the constant intercept in this limit, 
and it can therefore be distinguished from the impurity effects in an unique way. 



Experimental detection of the spin current 

In detecting the spin current by attaching a ferromagnetic electrode as in Fig. 3A, one 
should see a change of the electric current I depending on the direction of M. For n = 
10 19 cm~ 3 , the spin current density j* is estimated as 2eh~ 1 jy = (cr s /a)J z ~ J z . The spin 
current jy induces the electric current I when the magnetization of the electrode is along 
+ re-direction; no current is observed for — x-direction, when the ferromagnetic electrode 
is a half-metal. In realistic situation, the ratio of the electric currents I(+x)/I(—x) is 
determined by that of the tunneling probabilities for parallel and anti-parallel spins at 
the interface with the electrode. This ratio is still expected to be well larger than unity. 

Another method of experimental detection of the spin current is to measure circular 
polarization of light emitted via recombination with electrons as in Fig. 3B. With the 
steady electric current J z within the layer, the spin current injected along the y-direction 
will be recombined with the electrons supplied from the n-GaAs in the (In,Ga)As quantum 
well. This can be regarded as the p-n junction with the opposite voltage for opposite 
direction of the spin s x . Therefore the current along the circuit is different between up- 
and down-spins with different direction, and hence the finite current is observed for I. 

One can also detect the spin current by measuring the spins accumulated near the 
edges of the sample. The distribution of spin accumulation is determined by a balance 
between the spin current supply and the spin relaxation. A diffusion equation for hole 
spins can be written as 

ds x (y,t) ^ d 2 s x (y,t) dj*(y,t) s x {y,t) 

~dt D dy* ~ dy ^~ {blb) 

where D = [iksT /e is the diffusion constant. For simplicity, let us assume a semi-infinite 
sample extending for y < 0. Considering a steady state and the spin current being 
assumed to be jy(x) = jy0(—y), we can solve the above diffusion equation as 

sM=fy ] f^evk. (S16) 



This solution means that the spins are distributed within the length scale L = yDr s with 
the total amount jyT s per unit area. At room temperature, from r s = lOOfsec (S8) the 
length scale L is estimated as L = 4nm. We estimate the upper bound for the electric 
current density as 10 4 A/cm 2 , and the corresponding spin current 2ejy ~ 10 4 A/cm 2 . 
Therefore the spin density integrated over the depth is given by jyT s ~3x 10 9 /iBcm" 2 , 
too small to be observed optically. However there are several ways to overcome this 
difficulty. One is to lower the temperature. At T = 30K, r s reaches around t s ~ 30psec 
(SI 6), nearly 300 times larger than that of room temperature, and the total spin density 
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becomes ~ 10 12 /ZBcm~ 2 . Another possibility is to inject the spin into n-GaAs through the 
p-n junction. The spin lifetime of electrons in GaAs is lOOpsec, 10 3 times longer than that 
of holes, and by suppressing the D'yakonov-Perel' mechanism it could be even as long as 
2nsec (S17). Therefore the accumulation depth L becomes (0.1-l)/xm, comparable to the 
wavelength of the visible light, with the area density 10 12 -10 13 /zbchi _2 for the spins near 
the edges of the n-type GaAs. Therefore it can be detected by the Kerr rotation through 
surface reflection. 
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